Electrons scattering in the monolayer graphene with the short-range impurities 
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Scattering problem for electrons in monolayer graphene with short-range perturbations of the 
types "local chemical potential" and "local gap" has been solved. Zero gap and non-zero gap kinds 
of graphene are considered. The determined S-matrix can be used for calculation of such observables 
as conductance and optical absorption. 



During the last years much attention was payed to the problem of the electronic spectrum of graphene (see a review 
Two-dimensional structure of it and a presence of the cone points in the electronic spectrum make actual a 
comprehensive study of the external fields effect on the spectrum and other characteristics of the electronic states 
described by the Dirac equation in the 2-f-l space-time. We consider in this work the electrons scattering in the 2-1-1 
Dirac equation model of the monolayer graphene due to the short-range perturbations. We do not take into account 
the inter-valley transitions. Particular attention to this case stems from the effectiveness of short-range scatterers in 
contrast to the long-range ones: an effect of the latter is suppressed by the Klein paradox 0. Short-range potential 
impurities in graphene were considered in works 0] , [l] j • In our work 0] , a new model of the short-range impurities 
in graphene was considered taking into account the obvious fact that the Kohn-Luttinger matrix elements of the short- 
range perturbation calculated on the upper and lower band wave functions are not equal in a general case. This means 
that the perturbation must be generically described by a Hermitian matrix. We considered the diagonal matrix case 
corresponding to a presence of the potential and mass perturbation. The bound states dependence on the perturbation 
parameters was studied in [6;] within the framework of this model. 

In the present work we study the electrons scattering by the short-range impurities within the framework of the 
model suggested in Q. 

The Dirac equation describing electronic states in graphene reads 

^-ihvF 7^9^ ~ 7o ("1 + Sm) vlj ip = {E -V)ip, (1) 
where vf is the Fermi velocity of the band electrons, 7^ are the Dirac matrices 

70 = CTg, 7i = (71, 72 = 2(72, 

(7i are the Pauli matrices, 2towf^ = Eg is the electronic bandgap, t/j (r) is the two-component spinor. The electronic 
gap can appear in the graphene monatomic film lying on the substrate because of the sublattices mutual shift 
The spinor structure takes into account the two-sublattice structure of graphene. Sm (r) and ^(r) are the local 
perturbations of the mass (gap) and the chemical potential. A local mass perturbation can be induced by defects in 
a graphene film or in the substrate f^. We consider here the delta function model of the perturbation: 

(5m (r) = — M(r — tq), V(r) ^ —ad{r — ro), (2) 

where r and tq are respectively the polar coordinate radius and the perturbation radius. Such short-range perturbation 
was used in the (3-|-l)-Dirac problem for narrow-gap and zero-gap semiconductors in ISj]. The perturbation matrix 
elements 



diag{Vi,V2)S{r - ro) (3) 

are related to the a, b parameters as follows 

-Vi^a + b, -V2=a-b (4) 

The delta function perturbation is the simplest solvable short-range model. Finite radius tq plays a role of the 
regulator and is necessary in order to exclude deep states of the atomic energy scale. The finite perturbation radius 
ro leads to the quasi-momentum space form-factor proportional to the Bessel function that justifies our neglect of 
transitions between the Brillouin band points K and K' ^ . 
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Let us present the two-component spinor in the form 



exp {—iEt) 



' fj (r) exp[i(j - 1/2)^] 
9j (r) exp [i (j + 1/2) ip] 



(5) 



where j is the pseudospin quantum number; j = ±1/2, ±3/2, .... In opposite to the relativistic theory, this quantum 
number has nothing to do with the real spin and indicates a degeneracy in the biconic Dirac point. The upper fj (r) 
and lower gj (r) components of the spinor satisfy the equations set 



^ + ^-9, -{E-m) f, ^{a + b) 6{r ~ ro)/„ 



(6) 



- + ~{E + m) g, = [a - b) 5{r - ro)g,. 
ar r 

These equations have a symmetry: 

/j ^9], E ^ -E, j - j, a -a. 
Let us introduce the function ipj (r) = fj/gj- It satisfies the equation: 

l/[(a + 5)^2^(a-fe)] 
Integrating in the vicinity of r = tq 

^ 



+ S{r - ro) = 



lim 



Vjiro+S) 



"5^0 J^, (ro^S) {a + b)(pj + {a- b) 



(7) 



(8) 



(9) 



(10) 



we obtain the matching condition 



arctan 



[if- ^/ {a + b) / {a - b)^ - arctan [if+ ^/ {a + b) / {a - b)^ = - b^, 



2' > b\ 



(11) 



where Lp^ = ip^ (ro — 5) , ip1j = ipj (ro + (5) , (5 — > 0. Excluding the spinor component gj from the equation set Eq. 
(|6|), Eq. ([7]) in the domains < r < rp and r > tq, we obtain the second-order equation: 



h = 0. 



(12) 



This equation is related to the Bessel one. We assume E to be real and satisfying the inequality E"^ > m^. Then the 
general solution of Eq. (fT^ in the region < r < reads 

(13) 



where n = ^/ E"^ — m? is the principal value of the root; (z) and N^, (2) are respectively the Bessel and Neu- 
mann functions. The constant C2 vanishes in the domain < r < ro since the solution must be regular at the 
origin. Expressing the (jj-component from Eq. ([7]), we can write 



9j = \l §Tp^^'^i'^J+i/2 i^r) ■ 



Thus 



E + m Jj-1/2 (Kro) 
E-m Jj+1/2 (kto) ' 



(14) 



Then we can obtain from Eq. (jlip : 
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arctan 



and, therefore, 



'fit (i^ro) = 



J,-i/2 {>^ro) -{a~b)T (a, b) J,+i/2 (^^0) 



Jj+1/2 i^ro) + ia + b) J§^T (a, b) J,_i/2 (^^o) 



tan( \/a^— ) .^9 ,9 



where T (a,b) is given by the formula: 

T(a,5) = < 

On the other hand, an expression for ip^ {i^tq) can be written similarly to 1141 



(16) 



(17) 



+ , , /; rBT^i/]!_\/2('*ro) + 5,i7]i\/2(^ro) 

^ ^ ^ V^i.S/,(..o)...<\/.(..o)' ^^^^ 

where S'j (k) is a phase factor of the out-going wave, i. e. S-matrix element in the angular momentum representation. 
Substituting Eq. (fT5)) into Eq. (fTB)) . we obtain an explicit expression for Sj (E): 

jr(2) 

SAE) = -^, (19) 

where 



(a) 



T(a,5) 



E — m 
E + m 



(a - 6) Jj+i/2 (Kr-o) ('^^o) + J ^ (a + 6) Jj-1/2 (Kro) H]^1'i/2 ('^^o) 



E + m 



r(") 



E — m 



■ (20) 



Here a takes values 0, 1. Since iji^^ (z) = hI^^* (z) for real z, the scattering matrix is unitary everywhere on the 
continuum spectrum. Eq. p9p solves the electron scattering problem for the given potential. The denominator of 
Sj(E) is just the left-hand side of the characteristic equation derived in 6*1. Imaginary roots of it correspond to the 
real energy eigenstates (bound states) lying in thegap, which were studied in that paper. The characteristic equation 
reads 



(21) 



(kJ-o) ^^]+1/2 ('«''o) - J3 + I/2 {l^fo) ^]-l/2 ('«^0)) 

T {a, b) J (a - b) Jj+1/2 (nro) i?j+i/2 (uro) + 



E + m 

E — 771 



(a + b) (ura) H^Z\/2 ('^'^o) 



(22) 



Using the relations Hn"^ (z) = J„ + iNn, H'iC' = Jn — iNn, we can write S-matrix in the form: 

^ A,{E)+7B,{E) ^ B, {E) + lA, {E) 
^ ^ ' Aj [E) - iBj [E) Bj (E) - iAj (E) ' 



r(2) 



(23) 
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and, therefore, it can be presented in the standard form [lO| 

Sj (E) = exp [i2Sj (E)] 
where the scattering phase is given by the expression 



Sj [E] = arctan 



A, [E) 

B, {Ey 



(24) 



(25) 



Formulae (|23p . (|24p show that the scattering matrix Sj (E) is unitary on the continuum spectrum. The functions 
Aj (E) and Bj [E] are determined as follows 



A,{E)^-T{a,b) 



, , , I E + m 



E — m 



J:i-i/2 (i^ro) + {a-b) 



E — m 
E + m 



Jj+i/2 (i^ro) 



(26) 



B,iE) = r(a,6) 



[jj+1/2 (^ro) Nj_i/2 (Kro) - Jj-1/2 (^ro) iV,+i/2 (^ro)] 



(27) 



It is seen from (P5|) . (I?fl) that all (5^ (i?) vanish , when a and 6 tend to zero, i. e. in the absence of a perturbation. 
Using the Bessel functions expansion ^] 



J„ (x) ^ (l/nl) {x/2r , 

- (r(?i)/7r) (2/j;)" for7i>0, 
(2/7r)log(7^a;/2) for ?i = 



(28) 
(29) 



we conclude that for the low-energy scattering ktq « 1, 6j (E) is small as (nro)^''^'^^^^ except of j = ±1/2. Here 
log 7^ is the Eyler-Mascheroni constant. In the case of small radius and low energy E we can neglect all higher 
angular momentum partial waves taking into account only phases Sj for j = ±1/2 : 



tan(5i/2 (E) = 



= -T{a,b) 



i<^ + b)^^+{a-b)^^{.r,r. 



[(«ro/2) f log {lE^r,/2) - i (2/«;ro)] + T [a, b) {a + b) J|±^f log {lE^r,l2) + (a - b) ■ / E^Elll 



(30) 



tan (5_ 1/2 (^^) 



-T(a,5) 



(a - &) 1 


yi^^ + (« + ^)v/i^(-V2)^ 


{2/nr^) - (/tro/2) | log {-iE^r^/2) 


+ T(a,fe) 





r(a,5) 



E — m 
E + m 



(a - 6) TT , Kro — >0 



(31) 
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We see that the phase is proportional to nro in the long- wave limit as it is necessary [13, The scattering 
amplitude / (9) and transport cross-section T,tr can be expressed in terms of Sj (E) as follow [3] : 

= — ^ ^ [S,iE)-l]exp[zU- 1/2)6], (32) 
T.tr = 2/k ^ sin^ {5j+i - 5j) (33) 

i=±l/2,±3/2,.. 

Near the resonance states the Breit-Wigner form of the phase is valid : 

5j w 5^^^ + arctan 



' ' 2(Ef^-E 

where e''^^ and T j are respectively the position and width of the resonance level, S'^p is the slowly-varying potential 
scattering phase. 

The presented above formulae can be used in order to calculate the Boltzmann conductivity [ll| 
where the transport relaxation time equals 

l/Ttr = NiVF^tr- (35) 

Here Ni is the areal impurity density, Ep = vfkf- The above equations transform a dependence of the scattering 
data on the Fermi energy and impurity perturbation parameters a and b into the correspondent dependence of the 
Boltzmann conductivity. Thus characteristic features of the scattering data determine a behaviour of the electric 
conductivity. Proper numeric calculations will be presented elsewhere. 
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